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We model combined photoassociation and Feshbach resonances in a Bose-Einstein condensate, 
where the shared dissociation continuum allows for quantum interference in losses from the conden- 
sate, as well as a dispersive-like shift of resonance. A simple analytical model, based on the limit of 
weakly bound molecules, agrees well with numerical experiments that explicitly include dissociation 
to noncondensate modes. For a resonant laser and an off-resonant magnetic field, constructive inter- 
ference enables saturation of the photoassociation rate at user-friendly intensities, at a value set by 
the interparticle distance. This rate limit is larger for smaller condensate densities and, near the Fes- 
hbach resonance, approaches the rate limit for magnetoassociation alone. Also, we find agreement 
with the unitary limit-set by the condensate size-only for a limited range of near-resonant magnetic 
fields. Finally, for a resonant magnetic field and an off-resonant laser, magnetoassociation displays 
similar quantum interference and a dispersive-like shift. Unlike photoassociation, interference and 
the fieldshift in resonant magnetoassociation is tunable with both laser intensity and detuning. Also, 
the dispersive-like shift of the Feshbach resonance depends on the size of the Feshbach molecule, 
and is a signature of non-universal physics in a strongly interacting system. 

PACS numbers: Pacs number(s): 03.75.Nt, 05.30. Jp, 34.50.Rk 



I. INTRODUCTION 

Photoassociation occurs when an ultracold pair of 
atoms absorbs a photon and jumps from the two-atom 
continuum to a bound molecular state [H |2] . Similarly, 
magnetoassociation occurs when one atom in an ultra- 
cold pair spin flips in the presence of a magnetic field 
tuned nearby a Feshbach resonance [3] H] . The two theo- 
ries are formally identical, and are referred to in general 
as coherent association. If the atoms are quantum degen- 
erate, then coherent association can be highly efficient, 
due to large phase space density E] . 

Pioneering experiments [7J with Raman photoassocia- 
tion of a Bose-Einstein condensate demonstrated highly 
efficient atom-molecule conversion that was just on the 
verge [8] of coherence. Subsequent experiments [9j [10] 
focused on the strongly interacting regime and the asso- 
ciated rate limit on atom- molecule conversion, as well 
as the free-bound lightshift. The rate limit for con- 
verting atoms into molecules is set by either unitar- 
ity [IT] , i.e., the size of the condensate, or rogue dissoci- 
ation [51Q21EI], i.e., the interparticle distance, and the 
lighthshift is predictedly red and proportional to laser in- 
tensity [Mj [15] . Whereas the observed lightshift in Bose- 
condensed Na [3] and ultracold 7 Li [TO] was consistent 
with theory [TU [13] and previous non-degenerate exper- 
iments |16j . the rate limit could not be reached in Na 
and, although it was reached in 7 Li, the absence of Bose 
condensation precluded a direct comparison between the 
rogue and unitary limits. Rate limit aside, an optically- 
tuned scattering length [T4J [17] has been observed in an 
87 Rb [18 condensate, along with observation [19] of the 
two photon atom- molecule dark state [T51 |2"0"] . 



On the other hand, the Feshbach resonance has gar- 
nered exceptional interest. Early experiments focused on 
tuning the s-wave scattering length in Na [3J, followed 
by the formation of a Bose-Einstein condensate of oth- 
erwise incondensable 85 Rb [21 and 133 Cs [22 . Later ex- 
periments with 85 Rb demonstrated the "Bosenova" col- 
lapse of a Feshbach-resonant condensate [23 , a counter- 
intuitive decrease in condensate losses for an increase in 
time spent near resonance [2~H |2"5] . as well as coherent 
oscillations between the Bosenova burst and its remnant 
core [25] [2H1 [27] . Soon after, magnetoassociation led to 
the molecular condensate milestone |28] and, in combina- 
tion with laser transfer of the vibrationally-excited Fes- 
hbach population [55], the production of a quantum de- 
generate gas of absolute ground state molecules 30J . The 
very latest results report that the Feshbach resonance has 
enabled exothermic chemical reactions in a quantum de- 
generate gas of dipolar molecules 31]. 

Meanwhile, a combination of the two resonances has 
also been investigated. One of the first experiments with 
Feshbach-resonant scattering lengths in 85 Rb used a pho- 
toassociation laser as a probe [35] , thereby observing en- 
hancement of the conversion from atoms to molecules, 
and later experiments in 133 Cs observed supression [33] . 
Supression and enhancement were then observed together 
in 133 Cs [3J], an investigation that also included a suc- 
cessful numerical model based on photoassociation with 
a Feshbach-tunable scattering length. Most recently, the 
Feshbach resonance was observed to vary the photoasso- 
ciation rate constant by over four decades, and to anoma- 
lously shift the position of laser resonance to the blue end 
of the photoassociation spectrum for near-resonant mag- 
netic fields [351 [35] . Finally, the most recent theory pre- 
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diets strongly enhanced production of dipolar molecules 
in the electro- vibrational ground state |37j , as well as an 
Autler-Townes-like splitting of the combined resonance 
lineshape for a high intensty photoassociation laser [35]. 
These models are analogous to laser-induced autoioniza- 
tion [39], and laser-enhanced photoassociation [40] , 

The purpose of this Article is to develop the details of 
the matter optics theory of combined photoassociation 
and Feshbach resonances. For a resonant laser and an off- 
resonant magnetic field, our simple analytical model [36 
for the rate constant and continuum shift of laser res- 
onance agrees well with numerical experiments, with 
the main result indicating that losses are stronger for 
a more dilute gas. Next, we connect the present results 
to collisional models, illustrating how the rate constant 
and lightshift have a node at the magnetic field value 
where the Feshbach-resonant scattering length matches 
the semi-classical size of the photoassociation molecule. 
We also show that, for combined resonances, the rogue 
limit for atom-molecule conversion agrees with the uni- 
tary limit only over a limited magnetic field range near 
the Feshbach resonance. While these results are indepen- 
dent of the laser parameters, we find that, for a resonant 
magnetic field and an off-resonant laser, the quantum 
interference and fieldshift depend on both the laser in- 
tensity and detuning. Moreover, in this case the nodes 
in the rate constant and lightshift occur at difference in- 
tensity/detuning position, and the fieldshift depends on 
the size of the Feshbach molecule. 

This Article is outlined as follows. In Sec. [IT] we first 
introduce the few-level diagrams, Hamiltonian, and cor- 
responding nonlinear Schrodinger equation for the com- 
bined resonance system. In Sec. |II A[ we develop an an- 
alytical solution to the nonlinear Schrodinger equation 
based on the assumption of weakly bound molecules. For 
one or the other of the molecular states is tuned off res- 
onance, Sec. |IIB | derives the equations of motion for the 
atomic and molecular probabilities, as well as the atom- 
molecule coherence which, in turn, leads to the rate equa- 
tion for losses from the condensate. 



(a) 
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we re- 



in Sec. 

port an algebraic solution to the nonlinear Schrodinger 
equation that can be implemented numerically at little 
computational cost, and describe our procedure for nu- 
merical experiments. Sec. II D presents the parameters 



for 7 Li. Results are presented in Sec. Ill Resonant pho- 
toassociation is detailed in Sec. |III A[ including the com- 
parison with numerical experiments, the connection to 
collision models, and the comparsion to unitarity. Res- 
onant magnetoassociation is covered in Sec. |III B[ with 
a focus on the role of the size of the Feshbach molecule. 
Conclusions are given in Sec. 
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II. MODEL 

Consider N atoms that have Bose condensed into, say, 
the zero momentum plane-wave (Kk — 0) state |0). Pho- 
toassociation and the Feshbach resonance then couple 
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FIG. 1: (color online) Few-level scheme for a condensate 
tuned nearby a combined photoassociation and Feshbach res- 
onance, (a) Basic three-level scheme, where a photoassoci- 
ation and Feshbach resonance couple the atomic condensate 
|0) and molecular the condensates |1) and |2), respectively, 
(b) A quasicontinuum accounts for dissociation to nonconden- 
sate levels, (c) Eliminating the noncondensate levels leads to 
an effective F-system, where the virtual continuum couples 
the two molecular states, and where the detunings Si include 
the free-bound redshift. (d) When the system is far from one 
resonance, magnetic or laser, the off-resonant molecular state 
can also be eliminated, leaving an effective two-level system, 
where the detuning v includes an anomalous Stark-shift. 



atoms in the state |0) to diatomic molecules of zero mo- 
mentum in the states |1) and |2), respectively. As per 
Fig. [TJa), this is the T^-system familiar from few-level 
laser spectroscopy. Annihilation of an atom (molecule) 
of mass m (M = 2m) from the atomic (ith molecular) 
condensate is represented by the second-quantized op- 
erator ciq = a (bi). This theory is the simplest, and 
molecules dissociate only back to the level |0). To be 
more complete, molecular dissociation to noncondensate 
levels should also be included [Fig. |TJb)] . These lev- 
els must be considered because a condensate molecule 
need not dissociate back to the atomic condensate, but 
may just as well create a pair of atoms with equal-and- 
opposite momentum, since total momentum is conserved. 
So-called rogue [5J H2] , or unwanted [H] , dissociation to 
noncondensate modes therefore introduces the operators 
a ±k . 

The second quantized Hamiltionian for this system is 
H = Hq + Hi , where the condensate Hamiltonian is 



Ho 

n 



[Sib$bi - Ki(b\aa + 0+0+6)] , (la) 

i=l,2 



and the noncondensate Hamiltonian is 

k i k^O 

(lb) 
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Here the bare photo-association (Feshbach) detuning is 
5[ (($2), the atom-molecule coupling is n\ (k 2 ), the ki- 
netic energy of an atom pair of reduced mass m, r = to/2 
is fi.£k = H 2 k 2 /2m r , and the momentum dependence of 
the coupling to the noncondensate modes is given by /ik 
(/2k) with fio = 1. Spontaneous decay of the photoas- 
sociation molecule has been neglected for now, whereas 
spontaneous decay of the Feshbach molecule [42\ is ne- 
glected once and for all. 

To obtain mean-field equations, the Heisenberg equa- 
tion for a given operator, ihx — [x,H], is derived from 
the Hamiltonian ([I]), and all operators are subsequently 
declared c- numbers. In a minimalist model, x repre- 
sents either the atomic (ith molecular) operator eik 
or the anomalous density operator = aka-k, where 
Ak arises from rogue dissociation to noncondensate atom 
pairs of equal-and-opposite momentum. Converting sum- 
mations over k to integrals over frequency e introduces 
the characteristic frequency ui p — hp 2 ' 3 /2m r . The re- 
sulting mean-field theory is 

ia = —Qia*bi — fl 2 a*b 2 , (2a) 
ik = 5% - \n t a 2 - J dsVefi{e)A(e), (2b) 
iA{e) = eA(e)-r2i/i(e)6i-n 2 / 2 (e)l>2. (2c) 



These amplitudes are of unit order, = V iV/tj oc ^fp, 

and £i = £li/(4ir 2 ujp^ 2 ). Finally, the continuum shape is 
related to the size of the underlying molecule and, since 
there are two molecules, we allow for different length 
scales with the Gaussian /j = exp[— (e/ (3i) 2 /2], where 
Pi = h/(2m r L 2 ) with L 4 the semi-classical size of the 
i-th molecular state. Here we choose a Gaussian over a 
Lorentzian |36j to foster faster numerical convergence. 



A. Weakly-Bound Molecules 

Previously, it was shown that the shared dissociation 
continuum can act like a virtual state [35J HO], leading 
to an effective cross-coupling between the two molecu- 
lar states. However, this result is based on adiabatic 
elimination of the anomalous amplitude which, despite 
capturing the essential physics, leaves something to be 
desired. This subsection therefore develops the cross- 
molecular coupling based on simple Fourier analysis. 

Consider the equations of motion ^ in the limit = 
with £ 7^ and f2j/i(s) 7^ 0. Fourier transforming 
the resulting equations and subsequently eliminating the 
rogue amplitude A(oS) gives 

ik(oj) = [ai + Ei(w)]6i(w) + !n 3 M62M (3a) 
ib 2 {u) = $ + £ 3 (w)]62(w) + ifl3(w)6i(u;), (3b) 

where the frequency dependent shift and cross molecular 



coupling are 

n 3 (w) = 



8tt 2 ( 



3/2 



47T 2 ( 



3/2 



J de\fe 



de\fe 



UJ — £ ' 

uj — e 



(4a) 
(4b) 



In the limit of weakly bound molecules, coi = 
luTi^oK[£i(w)], 7i/2 = lim^o^PiM], and f2 3 = 
lim^^o ^[^3( w )]- Taking the inverse transform then 
yields 

ik = (<Ji - »7i/2)6i - |n 3 62 (5a) 
ib 2 = (S 2 - i l2 /2)b 2 - (5b) 

where the continuum-shifted (renormalized) detuning is 
now Si = S[ - a 0i . 

The equations of motion (|5| imply that, to lowest or- 
der, the non-condensate Hamiltonian ( lb I is 



■±n 3 (b\b2 + blh), 



(6) 



Hi 

h ~ 2 " 

and the full mean-field theory (Oj 0) is, to lowest order, 
given as 

ia = — VL\a*bi — Cl 2 a*b 2 , (7a) 

ik = <?i6i - \n ia 2 - \n 3 b 2 , (7b) 

ib 2 = S 2 b 2 - \n 2 a 2 - ifi 3 6i. (7c) 

As per Eqs. (|6| and the shared dissociation con- 
tinuum of Fig.TTrb) can be approximated by an effec- 
tive cross coupling between the molecules, illustrated in 
Fig- 0c). 

Now, any time a continuum is coupled to a bound 
state, the eigenfrequcncy of the bound state is neces- 
sarily shifted. This shift is to the red, is well under- 
stood [5J EH Efl ED E3 j and has already been included in 
the renormalized detuning, cr i = 5ipi(0)] < 0. However, 
the continuum here is coupled to a pair of bound states, 
and the effective cross-molecular coupling means that the 
Feshbach bound state will shift the eigenfrequency of the 
photoassociation resonance, and vice versa. Returning 
to the Fourier expressions for the cross-coupled molecu- 
lar amplitudes [Eqs. ([3])], and simplifying the transform 
of the derivatives, yields the lowest order result for the 
shift of the photoassociation and magnetoassociation res- 
onances, respectively, 



0"! = 001 



2 = 002 



5 2 + S t ± ^(5 2 -S 1 ) 2 + n 2 
6x + S 2 ± ^J{5i-5 2 ) 2 + nl 



(8a) 



(8b) 



Here the positive (negative) root gives the shift above 
(below) resonance. In contrast to the shift for a sin- 
gle bound state, the bound Feshbach (photoassociation) 
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molecular state, along with the shared dissociation con- 
tinuum, leads to a net shift of the photoassociation (mag- 
netoassociation) resonance that turns from red to blue 
f° r <$2 — r2§/|cr i | (Si > — Ol/lcroal)) where we have as- 
sumed (ooi/03) 2 <C 1. In particular, the redshift depends 
on the size of the underlying molecule, and therefore the 
position of the node in the cross-molecular shift does too. 



B. Off- Resonant Systems 

We now detail a simple analytical model for the cross- 
coupled system where one of the fields, magnetic or 
laser, is tuned far from resonance. In this case, the off- 
resonant molecular amplitude can be adiabatically elimi- 
nated from the mean-field theory, leading to the effective 
two-level system shown in Fig. |ljd) . For the ith (jth) 
molecular state tuned onto (off) resonance, ibi/Si « 
and the mean-field theory Q becomes 



1 a 



-Xa*bi 



ibi = (y - iT/2)bi - \xa 2 \ 



(9a) 
(9b) 



where the effective atom-molecule coupling, the Stark- 
shifted detuning, and the effective damping, respectively, 
are given as 



X 
v 

r 



Si 

r, 



£jSj 



(9c) 
(9d) 
(9e) 



where Cj = 0§/4|^| 2 . Note that this lighthsift shift 
agrees with Eqs. ^ in the off-resonant limit Si ^> Sj, Q3. 
Here the decay rate of the photoassociation (Feshbach) 
molecule is Ti = T s + 71 = 72), where T s is the 
spontaneous electronic decay rate of the photoassociation 
bound state, spontaneous decay of the Feshbach state 
has been ignored, and the dissociative decay rate is 7$ = 
9[Ei(0)]. 

We now follow Ref. [5], and derive a rate equation for 
the atomic and molecular probabilities, Pq = a* a and 



Pi — 2b* bi, as well as the atom- molecule "coherence" 
Coi = aab*. Making extensive use of the product rule, 
e.g., iCoi — iaab* + 2iaab*, we have 

iPo = x(Ok-C&). (10a) 
iP % = -iTPi-xiCoi-CZi), (10b) 
l C 0t = -{ v + l T/2)C Qi + \ X Po{Po-2Pi). (10c) 



Solving Eqs. (JlOj) in the reservoir approximation, P a ~ 1 
and Pi <C 1, and for an adiabatic coherence Coi ~ 0, the 
rate equation for losses from the atomic condensate is 



P = pKPj. 



where the rate constant for condensate losses is 



pK 



1 



x 2 r 



(Ha) 



(lib) 



2 ^ 2 + r 2 /4' 

In addition to the anomalous lightshift, we realize that, 
as a function of the off-resonant detuning, Sj, the reso- 
nant rate constant K is suppressed below the rate for the 
on-resonant field alone, x < Oj, essentially vanishes for 
X ~ 0, and is enhanced above the rate for the on-resonant 
field alone for \ > Oj. Borrowing intuition from quantum 
optics, this suppression and enhancement is due to quan- 
tum interference between direct resonant association and 
association via the off-resonant molecular state. 



C. Numerical Experiments 

We also develop an exact numerical solution of 
Eqs. which is valid even for dual resonance. Off- 
resonant fields lead to stiff equations, so the usual nu- 
merical routines require long run times. On the other 
hand, most stiff routines are based on matrix approaches 
with run times that scale quadratically with the number 
of quasicontinuum states, leading to another impasse. 

Hence, we write the equations of motion |2]) in the 
matrix form hp = Hip, where semi-classical Hamiltonian 
and wavefuntion are defined by 



Hip 



Si 





5-2 





£1 

52 


/eT/i(ei)de 
/elf 2 [ex)d£ 




























£N Q 



h 

b 2 
a 

A( £l ) 



.(12) 



The molecular amplitudes have been scaled by \/2 so that the molecular probability is Pi — \bi\ 2 , and the inte- 



5 



gration over the noncondensate modes has been broken 
into a simple summation. Here we see that H is charac- 
teristically sparse, i.e., only the first two rows, the first 
two columns, and the diagonal contain non-zero elements. 
Hence, the double-resonance system can be solved using 
an adapted Crank-Nicholson routine [43] , with run times 
that scale linearly with the number of states. In partic- 
ular, given the solution V>o = V'(^) a t time t, the solution 
ip = tp(t + dt) at time t + dt is given as 



l-i(dt/2)H , 
^ = 1 + i(dt/2)H ^ 



(13) 



Defining M± = 1 ± [dt/2)H, Eq. g can be written ip = 
M-M^ipo so that, for <j> the solution to M + <p — ipo, then 
the desired solution is simply ip = M-<f). Nevertheless, 
H = H(ip) means that coherent association is inherently 
nonlinear, and a predictor-corrector approach is therefore 
taken: H [ipo] determines the solution ip from the initial 
solution ipo, and H[ip], in turn, determines the sought 
after solution ip. 

Finally, these are numerical experiments in the truest 
sense of the words. First, the numerical rate constant 
for photoassociation is defined by pKr = 1, where r is 
the time required for the condensate fraction to fall to 
| a. | 2 = 1/e « 0.3678. Also, the numerical lightshift is 
determined by minimizing r with respect to the laser 
detuning. In other words, we set the laser detuning to 
an inital guess Si) 



guess 



r s , determine r; reset the laser 
detuning, update r; repeating until r stops decreasing 
and begins to increase, which defines the minimum. 



D. Explicit Parameters 

We now determine typical parameters for each reso- 
nance. Starting with photoassociation, we set i = 1 and 
£2 = in Eq. (lib I. In the low-intensity limit, sponta- 



neous electronic decay dominates dissociation, 71 <C r s , 
since 71 oc f2 2 oc 7, so that the effective loss rate is 



r ps r s . Equation (lib I for the rate constant then be- 



comes pKo — 20,2/Ts- Defining Kq = K Q I and fix — 
fii vT, the photoassociation Rabi coupling is determined 
by fli = yj pK' a T s /2. In the high intensity limit, 71 ^> T s 
and decay is dominated by dissociation, r « 71, so that 
the saturated rate constant is defined by pK^ = 2Qf/7i, 
so that ji/I = KqT s / Kao. Accordingly, the saturation 
intensity is defined by I — Ts/{ji/I), given explicitly 
as J == Koo/K'o- Also, the characteristic length scale for 
saturated photoassociation is l\ = mK 00 /{&'KK). Lastly, 
the width of the photodissociation continuum is obtained 
from the lightshift per unit intensity uqi = Ooi^> so that 
fh = h/{mL\) and L x = mK'^T s / ' (iir 2 ha' 01 ) . 

Turning to magnetoassociation, experiments usually 
measure the resonant scattering length defined by a rcs — 
a^As/(B — Bo), where is the width and Bo is 
the magnetic field position of the Feshbach resonance, 
and where ctbk is the zero-field scattering length. In 



the present model, the resonant scattering length is de- 
fined from Eqs. Efy for i = 2 and C\ = 0, combined 
with the limit 62 ^> 72, so that a ICS — —H,^/^- Here 
we have assumed the continuum redshift is included in 
the definition of Bo- The Feshbach Rabi coupling is 
then = -\/87r/9|abk| AsAp/m, and the detuning is 
62 = (B — B )A fl /h, where A M is the difference in mag- 
netic moments between the Feshbach molecule and the 
free-atom pair. 

Focusing on 7 Li [10 , typical photoassociation pa- 
rameters are then T s = 12 x 2-k MHz, K' = 7.9 x 
10- 10 (cm 3 /s)/(W/cm 2 ), = 2.2 x 10~ 8 cm 3 /s, Si = 
1.7 x 2ir MHz/(W/cm 2 ), and p = 4 x 10 12 cm" 3 . The 
atom-molecule Rabi coupling is then fiiyTo = 290 x 
2ir kHz, and the photoassociation saturation intensity is 
Io = 28 W/cm 2 . From the formal definition of the pho- 
todissociation rate, 271 = 7r£f2 lv /w7 with u>x = h/(m£f), 
the photodissociation length scale is £\ = A/ 1.11, where 
2ttX = 671 nm is the wavelength of the photoassociat- 
ing light. Lastly, the width of the dissociation contin- 
uum is /3i = 29 x 2ir MHz, so that Li = 133a is con- 
sistent with the Condon radius of the photoassociation 
molecule [10j[35] (ao is the Bohr radius). Note that the 
cleaning up of dimensionlcss constants of order unity have 
led to slightly different values from Ref. [56"] . 

In 7 Li [35], a Feshbach resonance is located at Bo = 
736 G, the product of the zero-field scattering length and 
the resonance width is |obk|As = 1.6 nm-G, and the 
difference in magnetic moments between the Feshbach 
molecule and the free-atom pair is A M = 2po (the Bohr 
magneton is po)- The Feshbach atom-molecule coupling 
is then = 127 x 2ir kHz. The sole unfixed parameter is 
the kinetic energy of the magnetodissociated pair, Hu>2 = 
q|/2r7i r . For a zero-temperature homogeneous system of 
point-like Feshbach molecules, the reasonable ansatz for 
the length scale is the interparticle distance, so that the 
uncertainty principle gives (72 = hp 1 ^ 3 . 



III. RESULTS 

This section presents our analytical and numerical re- 
sults for resonant lasers and off-resonant magnetic fields, 
and vice versa. Off-resonant detunings are restricted to 
below threshold (negative) values, since condesates are 
unstable for negative resonant scattering lengths. 



A. Resonant Photoassociation 

1. Comparison to Numerical Experiments 

The purpose of this section is a detailed comparison be- 
tween our previously reported analytical model |36j with 
numerical experiments. Away from the Feshbach reso- 
nance and on lightshifted resonance (y = 0), the rate 
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constant is 



pK = 



9\ 



9% 



p 1 l z L\ ) \8irujp52 



-i 2 



91 



8ttu) p 5 2 



9.x 



-i 2 



72 



(14) 



where rj ~ 1 is a numerical factor leftover in leftover 
from the integration of the Gaussian continuum shapes, 
and we have assumed point-like Feshbach molecules 
L 2 <C L\ (when need for numerics, we use L 2 = 40ao). 
For increasing laser intensity, the saturated rate constant 
is 



pK = 



i + V 



i/a L ( %™ p 5 2 



2 «l 

72 



and the enhanced saturation intensity is 



To 



1 (P 1/3 L) 



9 2 
9qi 



roi 

72 



(15a) 



(15b) 



Near the free-bound resonance (5 2 ~ 0), the rate constant 
saturates at the rate constant for free-bound photoasso- 
ciation alone, pK — 9 2 /T 2 , and the saturation intensity 
tends to zero. On destructive interference, photoassocia- 
tion ceases for a magnetic field position that depends on 
the size of the Feshbach molecule: 



B 



in A m cj p 



Finally, the anomalous lightshift is 

1 



&i = cm 



1 



>r 



p 1 /' i LiJ \8irujpS 2 



91 



(15c) 



(15d) 



which has its own node at the magnetic field position 

(15e) 



'I 



h9\ 



B = B - — -^p 1/3 i!. 

S7T A p UJ p 



Results are shown in Fig. [2j As a function of laser in- 
tensity, the analytical and numerical rate constants agree 
best for low density, suggesting that the model is most re- 
liable when rogue dissociation to noncondensate modes 
is dominant. As a function of magnetic field, and for 
fixed intensity, the agreement holds steady except near 
resonance, most likely because the adiabatic approxima- 
tion used to eliminate the Feshbach molecular amplitude 
breaks down. The agreement between the analytical and 
numerical results for the lightshift is also good. 



2. Connection to Collision Models 

We may also make the connection to existing mod- 
els Hi] of photoassociation near a Feshbach reso- 
nance. For an off-resonant field, the resonant interatomic 
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FIG. 2: (color online) Rate constant and lightshift for res- 
onant photoassociation and off-resonant magnetoassociation, 
as a function of laser intensity (a,b) and magnetic field (c,d). 
Also in panel (a), dependence of the rate constant on density, 
where the blue solid (red dashed) line is the analytical result 
and the blue open circles (red diamonds) are the numerical 
results for p = 1 x 10 12 cm" 3 (p = 4 x 10 12 cm" 3 ). 



scattering length is defined 4Trhpa ICS /m = —9\/25 2 . The 
rate constant ( 14 1 is then 



pK 



L 



(16) 



72 



Similarly, the saturated rate constant, the saturation in- 
tensity, and the lightshift are re-written as 



pK 

IJn 
I01 
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Results are shown in Fig. [3j Both analytically and 
numerically, the rate constant for resonant photoasso- 
ciation ceases when the Feshbach-resonant scattering 
length matches the size of the photoassociation molecule, 
a r ~ L\. The node in the lighthshift occurs similarly. 
However, this coincidence is due to the limit where the 
Feshbach state is much smaller than the excited state, 
L 2 <^ L\ \ if the two molecular sizes are comparable (but 
not equal), or the Feshbach state were much larger, then 
the rate constant and the lightshift vanish at different 
magnetic field locations. 
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FIG. 3: (color online) Rate constant and lightshift for reso- 
nant photoassociation and off-resonant magnetic fields, as a 
function of the resonant scattering length. Here the density 
is p — f x 10 12 cm -3 , the solid line is the analytical result, 
and the open circles are the numerical results. 



3. Comparison to Umtarity 

Here we explicitly compare the rogue and unitary 
limits. The unitary limit is generally set by the 
de Broglie wavelength which, for a zero-temperature 
Bose-Einstein condensate, is given by the size of the con- 
densate. In the Thomas- Fermi model, the size of the 
condensate is determined from the trapping frequency 
v = {^1,^2,^3} and the condensate chemical poten- 
tial 2/z/cD = \J (15-/Va r0S /ah o ) 2 , where a^ = yjh/mv is 
the harmonic oscillator length scale with v — ^/T^v^Mz. 
Finally, a rcs is the Feshbach-resonant s-wave scattering 
length given by 4Trhpa ICS /m — ~ f2§/(2<52)- Focusing on 
cigar-shaped traps, the condensate is tightly confined in 
the a;y-plane and loosely confined along the z-direction, 
i.e., v\ = v<x = v r and v r 3> V3. The unitarity- limited 
rate constant is then 



K„ = 



hA d 



(18a) 



where the de Broglie wave length is set by the Thomas- 
Fermi radius of the cigar 



Ad = 2i?TF — 



8hp 



(18b) 



In Fig.gJ we see that the rogue and unitary limits only 
agree reasonably for magnetic fields tuned within a small 
window near the Feshbach resonance, and the two limits 
diverge considerably as the system moves away from and 
onto the Feshbach resonance. It is no surprise that uni- 
tarity does not account for the node, and we attribute 
the discrepancy near resonance to a breakdown of both 
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FIG. 4: (color online) Unitarity vs. rogue limit for reso- 
nant photoassociation and near-resonant magnetoassociation, 
as a function of magnetic field. The laser intensity is set to 
I = 100 W/cm 2 to ensure saturation absent Feshbach en- 
hancement, (a) The rogue (unitary) rate constant is given 
by the solid (dashed) line, (b) Ratio of rogue to unitary rate 
constants. The red diamond is not a numerical result, but 
simply marks B = 732 G for comparison to Fig. |5|a). 



the Thomas-Fermi and rogue models as the system ap- 
proaches threshold, where Feshbach molecules should be 
explicitly included. Far below threshold, where photoas- 
sociation alone dominates, the rogue limit is more strict 
than unitarity, as per Ref [12]. 



B. Resonant Magnetoassociation 

We also consider a resonant magnetic field and a de- 
tuned laser, in which case i = 2 and j = 1 in Eqs. (lib I. 
For fixed laser detuning and Oi = OiVT, the intensity 
position of the node is I c = (—8T:Siujp/rifli)/(p 1 / 3 Li). 



This node exists in addition to any nodes in the solu- 
tion to the radial Schrodinger equation [3], which oc- 
cur for arbitrary intensity. The details of the Feshbach 
resonance are largely irrelevant, and I c depends instead 
on the laser detuning and the classical size of the pho- 
tassociation molecule. Specifically, for a 7 Li condensate 
with p = 10 12 cm -3 , as well as a red laser detuning 
of one linewidth, S\ = —T s , the cessation intensity is 
I c w 25 W/cm 2 , as per Fig. [5] (a). Perhaps surpris- 
ingly, the rate constant actually increases as the intensity 
goes to zero, but this is only because magnetoassociation, 
the faster process, begins to dominate. Similarly, a node 
arises for fixed laser intensity and variable detuning, as 
shown in Fig. [5jc). Given the connection to the colli- 
sion model of the previous section, the dependence of 
the magnetoassociation node on both laser detuning and 
intensity arises because these are two means for optically 
tuning the resonant scattering length [T4l PTTl ITS] . 
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FIG. 5: (color online) Rate constant and fieldshift as a func- 
tion of laser intensity (a,b) and detuning (c,d), where the 
density is p — 1 x 10 12 cm -3 . (a,b) The blue solid (green 
dashed, red dotted) line is for laser detuning — (Ji/Toi = 4 (2, 
1). (c,d) The blue solid (green dashed, red dotted) line is for 
laser intensity I — 20 (40, 80) W/cm 2 . Note that a negative 
(positive) fieldshift means that magnetic field position of the 
Feshbach resonance decreases (increases). 



Beyond quantum interference, the photoassociation 
molecule shifts the magnetic field position of the Fes- 
hbach resonance, (72 = 002 + £fii/&3> &s shown in 
Fig. |5pt>,c). Besides the dependence on both laser inten- 
sity and detuning, it appears that the nodes in the reso- 
nant magnetoassociation rate constant and fieldshift oc- 
curs at very different values of intensity/detuning, in con- 
trast to the case of resonant photoassociation discussed 
in the previous section. This discrepancy is due to the 
difference in the sizes of the photoassociation and Fes- 
hbach molecules: if the molecular sizes were equal, the 
nodes would vanishes similarly; if the Feshbach molecule 
were taken as larger, then the magnetoassocaition nodes 
would match, and the photoassociation nodes would mis- 
match. Finally, the fieldshift node vanishes entirely if 
the cross-molecular shift, 81C/H3, never equals the back- 
ground shift (T02, i-e., if the optically-resonant scattering 
length never equals the size of the Feshbach molecule, as 
per Fig. [BT 



IV. CONCLUSIONS 

A quantum matter optics model of combined photoas- 
sociation and Feshbach resonances in a Bose-Einstein 
condensate has been developed. We reported the de- 
tails of an analytical model based on the weakly-bound 
molecule assumption, which is ultimately the same as 
assuming that the anomalous amplitude evolves adiabat- 
ically [T3J as well as exact numerical model that is 
valid even for dual resonance. For resonant lasers and 



off-resonant magnetic fields, the analytical results com- 
pare well with numerical experiments that account ex- 
plicitly for dissociation to noncondensate modes. This 
agreement is best for low density, suggesting that the 
weak-binding/adiabatic approximation is best applied 
to this regime, or that the magnetodissociation ansatz, 
72 oc /0 1 / 3 , breaks down. The combined-resonance rate 
constant is larger for smaller density, and the agreement 
between the analytical and numerical results indicates 
that this density dependence is inherent to the system, 
and not due to the magnetodissociation ansatz. 

Despite the agreement with numerical experiments on 
Feshbach-enhanced photoassociation, both approaches 
underestimate the saturated rate constant, K x = 1.4 x 
1CP 7 cm 3 /s, that is observed experimentally for resonant 
photoassociation at B = 732 G [35]. This discrepancy 
has two likely causes. First, the position of the node 
in the rate constant depends on the semi-classical size 
of the photoassociation molecule, L\ = 133cto, and we 
have over-estimated L\ by about 25% compared to the 
observed value, L x = 106ao [35] . A smaller value of 
Li would mean that B — 732 G is further from the 
node, and that enhancement is therefore stronger, mak- 
ing for a larger rate constant. Second, the combined rate 
constant increases for decreasing density, which means 
that the edges of the condensate, rather than the center, 
contribute more heartily to molecule formation in this 
regime. Better modeling of the molecular size, as well as 
full accounting of condensate inhomogeneity, could im- 
prove the agreement with observation. 

Considering the rate limit for Feshbach-enhanced pho- 
toassociation, we find that the rogue and unitary limits 
agree only in a limited range near the Feshbach reso- 
nance. This is not terribly surprising, since unitarity does 
not account for quantum interference, and both models 
breakdown on-resonance. For B = 732 G the two limits 
differ by about a factor of two, and nearer to resonance 
resonance the two results agree reasonably. This partic- 
ular disagreement is understood as a combination of an 
inaccurate size for the photoassociation molecule and the 
assumption of a homogeneous condensate. Moreover, it 
makes sense that a unitary limit based on a Feshbach- 
resonant scattering length agrees with our results near 
B-threshold where magnetoassociation losses dominate. 
Lastly, our previous [H] prediction for the photoassocia- 
tion rate limit is ~ 6 ui p , which compares favorably with 
both the 732 G result ~ 8w p , and the Feshbach-resonant 
rate ~ 9w p . Nevertheless, this agreement is to be taken 
with a grain of salt: the combined resonance rate limit 
saturates at the rate for magnetoassociation alone (times 
an interference factor), which need not be the same as 
the rate limit for photoassociation alone. 

Finally, we discuss the result for resonant i?-fields and 
off resonant lasers. In particular, we focus on the impli- 
cations for many-body universality, whereby the only rel- 
evant length scale in a strongly-interacting system is the 
interparticle distance. In photoassociation, the existence 
and size of the bound molecule state is well established, 
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even in a strongly interacting (high laser intensity) sys- 
tem, and its appearance in Feshbach- assisted photoasso- 
ciation is therefore not terribly surprising. On the other 
hand, because many-body concepts dominate Feshbach 
resonance models, so too does universality, and the Fes- 
hbach molecule is generally neglected. Here we see that 
a node in the net shift of the laser-assisted magnetoas- 
sociation resonance is an important test of the validity 
of universal models of the Feshbach resonance. In par- 
ticular, the zero-field shift is a i oc 1/L 2 , which tends to 
infinity for point-like Feshbach molecules, and the node 
in the ficldshift then ceases to exist. A finite node in the 



fieldshift therefore indicates a lengthscale other than the 
interparticle distance, and a breakdown of universality. 
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